Thermodynamic properties of nuclear matter with three-body forces 
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We calculate thermodynamic quantities in symmetric nuclear matter within the self-consistent 
Green's functions method including three-body forces. The thermodynamic potential is computed 
directly from a diagrammatic expansion, implemented with the CD-Bonn and Nijmegen nucleon- 
nucleon potentials and the Urbana three-body forces. We present results for entropy and pressure 
up to temperatures of 20MeV and densities of 0.32 fm~'^. While the pressure is sensitive to the 
inclusion of three-body forces, the entropy is not. The unstable spinodal region is identified and the 
critical temperature associated to the liquid-gas phase transition is determined. When three-body 
forces are added we find a strong reduction of the critical temperature, obtaining Tc — 12 MeV. 

PACS numbers: 21.30.Fe, 21.65.Mn 



The determination of thermodynamic properties of hot 
nuclear matter is extremely relevant for its applications 
in astrophysics and heavy-ion reactions. Pressure and 
entropy at finite temperature are crucial ingredients in 
the modeling of core-collapse supernovae and protoneu- 
tron stars, which involve processes characterized by den- 
sities up to several time the nuclear saturation density 
and temperatures up to few tens of McV's. The nuclear 
equation of state (EOS) plays an important role also in 
the interpretation of nucleus-nucleus collisions, in which 
a hot and dense state of matter is formed. Of partic- 
ular interest is the study of a possible phase transition 
occurring at subsaturation densities. Because of the van 
der Waals nature of nucleon-nucleon (NN) interactions, 
at low densities and high temperatures nuclear matter 
is expected to undergo a first-order transition to a gas 
phase [ll]. It is claimed that an evidence for such behav- 
ior is found in intermediate energy heavy ion collisions, 
as a plateau of the caloric curve derived for the light 
fragments as a function of the reaction energy 0]. 

The liquid-gas transition in nuclear matter presents 
some differences with respect to the case of finite nuclei 
where, as the temperature increases, Coulomb forces to- 
gether with the decrease of the surface tension trigger 
the onset of mechanical instabilities. One can try to re- 
late the limiting temperature of nuclei T/ to the critical 
temperature in infinite matter Tc by properly taking into 
account Coulomb and surface effects [3|. Different esti- 
mates from effective models based on Skyrme forces Q 
and more microscopic approaches Q are in agreement 
and yield a ratio T//Tc of about 1/3. 

A reliable many-body theory is called for the interpre- 
tation of the heavy ion reactions in the goal of extracting 
the equation of state and the critical point of the nuclear 
liquid-gas phase transition Q and in order to extrapo- 
late to higher densities and arbitrary isospin asymmc- 
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tries. There exist only few realistic calculations of the 
nuclear matter EOS at finite temperature. The varia- 
tional approach, which yields reliable results at T = 0, 
is usually extended to finite temperatures by neglecting 
the modifications of the correlation functions Q . There 
are attempts of developing a finite T variational tech- 
nique but no calculations of thermodynamic quan- 
tities are available so far. The Briickner-Hartree-Fock 
method generalized by means of the Bloch-dc Domini- 
cis formalism has been applied to finite temperature nu- 
clear mat terjiol. [ll| . The Green's functions in-medium 
T- matrix [12, ^M, approach is suitable for computing con- 
sistently microscopic properties and thermodynamic ob- 
servables and naturally takes into account finite temper- 
ature correlations [l^ . Thermodynamic relations such as 
the Hugenhotlz-van Hove and Luttinger identities [lE,|TB| 
are automatically fulfilled by the ^-derivable T-matrix 
approximation (itI. [Tsj. 

In two recent papers we presented the first calculations 
of entropy and pressure at finite temperature within 
the thermodynamically consistent T-matrix approach 
and the first results which include three-body forces 
in the finite temperature Green's functions method in 
nuclear matter [20|. Following these works we compute 
here pressure and entropy without and with three-body 
forces in the case of symmetric nuclear matter for two 
different realistic NN potentials, in particular addressing 
the unstable region related to the liquid-gas phase 
transition. 



Let us briefly recall the adopted computational 
method, whose details can be found in Refs. [ll,[2|. We 
solve the set of coupled equations which involve the cal- 
culation of the T-matrix, summing the ladder diagrams 
at all orders, the determination of the nuclcon self-energy 
and the single-particle Green's function. The full off-shell 
propagation of particles in the medium is taken into ac- 
count. We employ two different realistic nucleon-nucleon 
interactions, the CD-Bonn [2l| and the Nijmegen [23| pa- 
rameterizations, together with the semi-microscopic Ur- 
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bana three- nucleon potential [23|. Three-body forces, 
necessary for a rehable description of the saturation prop- 
erties in symmetric nuclear matter, are included via an 
effective two-body interaction derived after averaging out 
the third particle. The effect on the two remaining nucle- 
ons results in a mean field whose two parameters are fixed 
by requiring the correct saturation density and binding 
energy. The range of temperatures and densities stud- 
ied is well above the superfluid transition in symmetric 
nuclear matter (23 |. 

After calculating the single-particle propagator G and 
the nucleon self-energy E from the iterative scheme of the 
T-matrix equations one can compute the grand canonical 
potential, expressed as 



n = -Tr{ln[G"^]} - Tr{SG} + $ . 



(1) 



The generating functional $ has a diagrammatic expan- 
sion similar to the one for the interaction energy in the 
T-matrix approximation, differing by a factor \/n [n is 
the number of interaction lines) in front of each diagram 
(see Fig. [T]). Therefore we derive $ from the expression 
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FIG. 1: Diagrammatic expansion of the generating functional 
$. The sum includes ladder- type diagrams up to infinite order 
(n is the number of interaction lines in each diagram). 



for (Hint) by multiplying the potential y by a numerical 
factor A and integrating 



1 dX 

T 



(2) 



where G is the dressed single-particle propagator com- 
puted in the system with the full strength interaction. 
The expectation value {Hint) itself can be easily obtained 
from the T matrix [l9| . From the grand canonical poten- 
tial (dividing by the volume V) we obtain the pressure 
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which is then computed directly from its diagrammatic 
expansion without the use of numerical derivatives. 

In Figs. [5] and [3] we present the pressure for symmet- 
ric nuclear matter without and with three-body forces 
for the CD-Bonn and for the Nijmegen potential respec- 
tively. Results are shown for different temperatures as 
function of the density. When looking at the curves ob- 
tained with the two-body NN potentials (upper panels) 
we notice that they cross the thermodynamically forbid- 
den region in which the derivative of the pressure with 
respect to the density is negative. For both CD-Bonn and 
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FIG. 2: (Color online) Pressure as function of density in sym- 
metric nuclear matter for different temperatures with the CD- 
Bonn potential (upper panel) and the CD-Bonn potential plus 
three-body forces (lower panel). 



Nijmegen interactions the inclusion of three-body forces 
(lower panels) leads to a stiffening of the dependence of 
the pressure on density. The spinodal instability region 
below saturation density is still present but gets substan- 
tially reduced. 

In order to address more in details the critical behavior 
and the effect of three-body forces on the critical tem- 
perature for the liquid-gas phase transition we study the 
limits of the instability region in the density-temperature 
plane. The spinodal region is characterized by the nega- 
tive derivatives of the pressure and the chemical potential 



dP 
d p 



<0, 



d n 
dp 



< 



(4) 



3 



T = 8 MeV 
T= 10 MeV 
T= 12 MeV 
T= 14 MeV 
T= 16 MeV 
T= 18 MeV 
T = 20 MeV 




-1 - 



I ... I ... I ... I ... I ... I ... I ... I ... I 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 
P/Po 



FIG. 3: (Color online) Pressure as function of density in sym- 
metric nuclear matter for different temperatures with the Ni- 
jmegen potential (upper panel) and the Nijmegen potential 
plus three-body forces (lower panel). 



Inside this region the system is unstable and tends to 
separate in two different phases, gas at a lower and liquid 
at a higher density. The two phases may coexist in an 
interval of densities and temperatures up to the point in 
which 
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FIG. 4: (Color online) Coexistence and spinodal lines for the 
CD-Bonn potential without and with three-body forces. 
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FIG. 5: (Color online) Coexistence and spinodal lines for the 
Nijmegen potential without and with three-body forces. 



P{Pgas) = ^'(Pliquid) (5) 

and 

M(Pgas) = M(Pliquid) • (6) 

The two regions end up in a coincident point which de- 
fines the critical density and temperature of the liquid- 
gas phase transition. 



We compute the spinodal and the coexistence lines for 
the CD-Bonn and the Nijmegen potential without and 
with the inclusion of three-body forces. The limits of 
the spinodal region defined by Eq. ^ are found as local 
maxima and minima of the pressure. The conditions on 
the pressure and chemical potential are equivalent for a 
one-component system such as symmetric or pure neu- 
tron matter. For the case of arbitrary isospin asymme- 
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tries one should in principle consider both of them, but 
there is evidence that chemical and mechanical instabil- 
ities coincide also for asymmetric nuclear matter psj . 

We use two methods for determining the spinodal line: 
we compute the pressure derivative directly from the in- 
terpolation of P and apply the one suggested by Baldo 
and Ferreira which consists in plotting the chemical 
potential as a function of the pressure, looking for the 
back bending of the curves. The results obtained with 
the two techniques do coincide, confirming the stability 
of the calculations even in the low density regime. The 
critical lines are displayed in Figs. [Hand [5] for CD-Bonn 
and Nijmegen respectively. At low temperatures the un- 
stable and coexistence phases are present over a large 
range of densities. As the temperature increases the in- 
stability region becomes smaller ending up with the point 
defining the critical temperature for the liquid-gas tran- 
sition. The critical temperature is Tc = 18 MeV for the 
CD-Bonn and = 20.5 MeV for the Nijmegen poten- 
tial. Both the spinodal and the coexistence regions get 
strongly reduced when three-body forces are included, 
with the critical point at Tc = 12.5 MeV when CD-Bonn 
is employed and Tc — 11.5 MeV if we consider the Ni- 
jmegen potential. 

Another quantity which characterizes the liquid-gas 
phase transition is a dimcnsionless parameter that com- 
prises the critical pressure, density and temperature, 
Pc/{pcTc), which assumes the value 3/8 for a van der 
Waals equation of state. In Table [J we summarize these 
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Tc (MeV) 


Pc (fm-3) 


Pc (MeVfm-3) 


Pc 

pc Tc 


CD-Bonn 


18 


0.107 


0.43 


0.22 


CD-Bonn + TBF 


12.5 


0.096 


0.14 


0.12 


Nijmegen 


20.5 


0.094 


0.50 


0.26 


Nijmegen + TBF 


11.5 


0.088 


0.15 


0.14 



TABLE I: Critical values of the temperature, density and 
pressure of the liquid-gas phase transition for the CD-Bonn 
and Nijmegen potentials without and with three-body forces. 



the Bloch-de Dominicis finite temperature generalization 
of the Briickner-Hartree-Fock method. Using the Ar- 
gonne VIS potential they estimate the critical temper- 
ature to be Tc ~ 21 MeV. However, when they include 
the Urbana three-nucleon potential in the calculations 
they find, in contrast with our result, that three-body 
forces do not strongly affect the critical temperature, 
which is reduced to Tc ~ 20 MeV. Other calculations 
within the Bloch-de Dominicis formalism [2^ or the rel- 
ativistic Dirac-Briickner-Hartree-Fock approach 0, [29l |. 
on the other hand, yield lower values of Tc, respectively 
9, 12 and 10.4 MeV, closer to our result with three-body 
forces. 

We compute the entropy per particle in the interacting 
system from the thermodynamic relation 
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The results for symmetric nuclear matter at p 
are shown in Fig. [S] for different temperatures. 
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critical quantities for the various potentials without and 
with three-body forces. The values for the CD-Bonn and 
the Nijmegen potentials are rather similar, with a strong 
decrease of the dimcnsionless parameter after the inclu- 
sion of three-body forces, signaling a departure from the 
van der Waals equation of state. 

Our estimates with the two-body nuclcon-nuclcon in- 
teractions only are in agreement with other microscopic 
calculations. Rios ct al. [2^ use a self consistent Green's 
functions approach but compute the free energy from 
the Carneiro-Pcthick quasiparticlc entropy p7j . obtain- 
ing for the CD-Bonn potential a critical temperature 
Tc = 18.5 MeV and Pc/{pcTc) = 0.20. When using the 
Argonne VIS parameterization, however, they find that 
Tc is reduced to about 11 MeV and Pc/{pcTc) = 0.14. 
Baldo and Ferreira performed calculations within 



FIG. 6: (Color online) Entropy per particle for symmetric 
nuclear matter at p = po with the CD-Bonn and Nijmegen 
potentials without and with three-body forces. 



effects of three-body forces are very small, as well as 
the dependence on the nucleon-nuclcon potential. The 
entropy appears to be independent of the details of the 
interaction, supporting the conclusion that it is not 
much affected by nucleon correlations. The reliability of 
calculations of the entropy at the level of a quasiparticle 
approximation has been checked also by Rios et al. 
[30[ using the Carneiro-Pethick dynamical quasiparticle 
formula 1271. 
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Existing estimates of the critical temperature for the 
hquid-gas transition in infinite nuclear matter are in the 
range Tc ~ 9 — 21 MeV. Within the finite temperature 
Green's functions approach we calculate thermodynamic 
properties of symmetric matter and study the spinodal 
instability region related to the first-order phase transi- 
tion. We observe that three-body forces have a strong 
effect on the pressure of the interacting system, reducing 
the size of the unstable region in the density-temperature 
plane. With the CD-Bonn and the Nijmegen potential 
we find that the critical temperature decreases respec- 
tively from Tc = 18 MeV to = 12.5 MeV and from 
Tc = 20.5 MeV to = 11.5 MeV. In all cases the 
critical density lies in the range pc — 0.09 — O.llfm"'^. 
The pressure at which the phase transition takes place 
decreases from Pc ~ 0.43 — O.SOMcVfm^'^ to about 
Pc = 0.15 MeV fm^^ after the inclusion of the threc- 



nucleon potential. 

Three-body forces do not affect the entropy, which 
turns out to be independent of the potential used in the 
calculations. This last results confirms that entropy is 
not much sensitive to nucleon-nucleon correlations (cf. 
Refs. [illS^). 
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